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Abstract 

We add a periodic potential to the simplest gravitational model of a superconductor and 
compute the optical conductivity. In addition to a superfluid component, we find a normal 
component that has Drude behavior at low frequency followed by a power law fall-off. Both 
the exponent and coefficient of the power law are temperature independent and agree with 
earlier results computed above Tg. These results are in striking agreement with measurements 
on some cuprates. We also find a gap A = 4.0 Tc, a rapidly decreasing scattering rate, and 
"missing spectral weight" at low frequency, all of which also agree with experiments. 
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1 Introduction 

Gauge/gravity duality has provided a powerful new tool to analyze certain condensed matter sys- 
tems. In particular, it can be used to calculate transport properties of strongly correlated systems 
at finite temperature. Remarkably, this is achieved by mapping the problem to a gravitational 
problem in one higher dimension. Since the system of interest lives in one dimension less than the 
gravitational problem being solved, this approach is often called holographic. 

The optical conductivity in a simple holographic model of a 2 + 1 dimensional conductor was 
recently studied including the effects of a lattice [1]. Earlier studies assumed translation invariance 
which implied momentum conservation. In that case, the charged particles cannot dissipate their 
momentum so the real part of the optical conductivity always contains a delta function at zero 
frequency reflecting infinite DC conductivity. With the lattice included, the delta function is 
resolved. It was found that at low frequency the conductivity follows the simple Drude form, but 
at intermediate frequency, it follows a power law 

kMI = ^ + c (1.1) 

The exponent —2/3 is surprisingly robust, and is independent of the all the parameters in the 
model including lattice spacing, lattice amplitude, and temperature. Strikingly, a power-law at 
mid-infrared frequencies with exactly this exponent is seen in measurements of the normal phase 
of bismuth-based cuprates [21 [3]. (See also |1] for earlier experimental results and [5], [Gj [7] for 
attempts to explain this scaling). Also in agreement with the measurements is the fact that the 
coefficient B is temperature independent, although the cuprates do not appear to have the off-set 
C. 



The origin of (1.1) from the gravitational side is still not understood. It has been shown that a 



similar power-law is seen in the thermoelectric conductivity and for a 3 -|- 1 dimensional conductor 
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(although the exponents differ) [8]. In a recent paper, Vegh [9] studied the optical conductivity 
in a holographic model in which momentum conservation was broken by a graviton mass term. 
He finds a similar power-law scaling of the conductivity (including the constant off-set), but the 
exponent is not fixed. It depends on the mass of the graviton. 

The goal of this paper is to extend these results to the superconducting regime. We add a 
lattice in the form of a periodic potential to the simplest (original) holographic superconductor 
[T0| [TT] and compute the optical conductivity. We find that below Tc, in addition to the superfiuid, 
there is a normal fluid component with the following properties 

• It has simple Drude behavior at low frequency 

where the normal component density p„ and relaxation time r are temperature dependent 
(but frequency independent). 

• The scattering rate 1/r drops rapidly below T^.. 



• At intermediate frequency, the normal component again satisfies the power law (1.1) with 
the same coefficient B that was seen above T^.. 

• There is evidence for a superconducting gap of size A = 4.0 T^. 

• Despite this, in the limit of low temperature, p„ does not vanish, indicating the presence of 
uncondensed spectral weight. 

• The Ferrel-Grover-Tinkham sum rule relating the superfiuid density to the decrease in 
Re[(T(a;)] is satisfied, but only if large frequencies of order the chemical potential are included. 
If one considers only the Drude peak and power law region, there is "missing spectral weight" . 

As we will discuss, these are all observed properties of the bismuth-based cuprates. In particu- 
lar, experiments have shown that the power-law is unchanged when the temperature drops below 
the critical temperature and the material becomes superconducting [12]. 



2 Gravitational model 



We will work with the simplest holographic superconductor which requires gravity coupled to a 
Maxwell field and charged scalar $. The action is 
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(2.1) 



where L is the AdS length scale and F = dA. The scalar mass, = —2/L'^, is chosen since for 
this choice, the asymptotic behavior of $ is simple. From here on we work in units in which L = 1. 
If the metric asymptotically takes the usual form 

-dt^ + dx^ + dy'^ + dz"^ 



ds' 



(2.2) 
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then 

$ = ^01 + 2^02 + 0{z^) . (2.3) 

In the standard interpretation, $ is dual to a dimension two charged operator O in the dual theory 
with source 0i and expectation value 02- Since we want this condensate to turn on without being 
sourced, we set 0i = 0. 

As Gubser first suggested p!3], the above action has the property that (electrically) charged 
black holes become unstable at low temperatures to developing scalar hair. The reason is essentially 
that the effective mass of the scalar gets a contribution e^A^g^^ < from its coupling the Maxwell 
field which causes the $ = solution to become unstable. It was shown in [TOl [TT] that this is 
precisely the gravitational dual of a conductor /superconductor phase transition. 

The vector potential At must vanish at the horizon and is asymptotically 

At = ii-pz + 0{z^) (2.4) 

where is the chemical potential and p{x) is the charge density. For these electrically charged 
solutions, the field equations and boundary conditions require the phase of $ to be constant. We 
will set it to zero and treat $ as a real field. We introduce the lattice by requiring that the chemical 
potential be a periodic function of x: 

fi{x) = + Aq cos{kox)] . (2.5) 

The lattice is only introduced in one direction for computational convenience, and we will compute 
the conductivity only in the direction of the lattice. This "ionic" lattice has been discussed earlier 
(see, e.g., [m [151 [161 [TTl HE] ) but almost always treated perturbatively. In [8] it was treated exactly 



in the theory (2.1) without the charged scalar field. 

Our model has several parameters. In principle one could vary the mass and charge of the scalar 
field $ in the action. However, we have already fixed the mass for convenience. The boundary 



condition (2.5) depends on three parameters: the mean chemical potential /i, the lattice wavenum- 
ber ko and the lattice amplitude Aq. In addition, the solution will depend on a temperature T. 
Due to a scaling symmetry, physics depends on only three dimensionless quantities which can be 
taken to be Aq, ko/p,, and T//2. 

For definiteness, we will choose ko/fl = 2 for the results presented later in this paper. Although 
the physics is scale invariant, when doing calculations, we will also set ft = 1. Our choice of ko 
is very close to the one which yields a low temperature DC resistivity which is linear in T in the 
absence of superconductivity [I11[I]. Even though feg corresponds to an unrealistically large lattice 
spacing, the nonlinearity of Einstein's equation generates structure on much smaller scales. 
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3 Backgrounds 



To numerically construct the gravitational dual of the normal phase we employ the ansatz used in 



-Hi G{z) (1 - z)dt^ + . + Si {dx + F dzf + ^2 dy'' 



A = i)dt, $ = 0, 
(3.1) 

and G{z) = 1 + z + z"^ — ^\z^ /2. Here G{z) controls the black hole temperature given by T = 
G(l)/47r = (6 — fil)/8n and H12, Si^2, F and ip are six functions of x and z, that we determine 
using the numerical methods described in P, which were first introduced in [I9] and studied in 
great detail in [20] . 

When = the problem reduces to the original holographic superconductor [TT]. It was 
shown there that there is a critical value of T//i such that above this value, the scalar field 
vanishes and the solution is the simple planar Reissner-Nordstrom AdS metric. Below this value, 
the Reissner-Nordstrom solution becomes unstable to developing scalar hair. 

We now increase Aq and solve the Einstein-Maxwell equations to find rippled versions of the 
Reissner-Nordstrom AdS metric. This is exactly the same as what was done for the ionic lattice 
in [8]. To see when these solutions become unstable to forming scalar hair, we look for a static 
normalizable mode of the scalar field. It is intuitively clear that if one increases the charge e on 
the scalar field, it becomes easier for this field to condense and the critical temperature becomes 
higher. To find the critical temperature for a given charge it is convenient to turn the problem 
around. For each Aq and any temperature T, we find the value of the charge e such that T is 
the critical temperature for a field with charge e. At the onset of the instability, $ is a static 
normalizable mode of the charged scalar field and can be treated perturbatively: 

{V- -ieA''){Va-ieAa)^ + 2^ = (3.2) 

with $ being a function of x and z only. Also, the connections V and A are evaluated on the 
rippled Reissner-Nordstrom AdS black holes. Finally, one can recast the former equation to take 
the following appealing form 

_ V'^Va* - 2 $ = (-A^A'^) $, (3.3) 
which one recognizes as a positive selfadjoint eigenvalue problem for e^. At the boundary we 



demand $ to decay as in Eq. (2.3) and at the horizon we demand regularity. 

The results are shown in Fig. [Tj There are seven curves on this plot denoting seven different 
values of the lattice amplitude Aq between and 2.4. Each curve shows the expected rise in critical 
charge with temperature (or critical temperature with charge). Comparing the different curves for 
a given charge, we see that increasing the lattice amplitude also increases the critical temperature. 
This was first noticed in [5T] when the lattice was treated perturbatively. It can be understood as 
follows. As we said, the critical temperature of a uniform holographic superconductor is propor- 
tional to /i. By making vary periodically, one raises the maximum value of /i which induces the 
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Figure 1: For each T//i we plot the charge of the scalar field for which T would be the critical 
temperature. This is repeated for several values of the amplitude of the lattice. From the top 
down the lines represent: Aq = 0.0, 0.4, 0.8, 1.2, 1.6, 2.0, 2.4. Setting e = 2 we read off the critical 
temperatures used in this paper. 



scalar field to condense at a higher temperature. It then leaks into the regions where /i is smaller. 
This effect was seen in models of a holographic Josephson junction |22]- We will set e = 2 since 
this is a natural value for a superconducting condensate. 

Having found the critical temperature, one can next solve the coupled Einstein-Maxwell-scalar 
equations to find the solutions for T < Tc. In order to find the rippled superconducting phase we 
have to change the ansatz (3.1 ). The reason being that Eq. (3.1 ) is adaptec|^to probe geometries for 



which the entropy is nonzero as T — )■ 0. The homogeneous holographic superconductors are known 
not to have this feature. In particular, their entropy decreases towards zero as the temperature 
is lowered, thus obeying to the third law of thermodynamics. One ansatz that is adapted to this 
property is given by 

ds^ = — l -Hi y\ (1 - z)dt^ + , ' , + y\ \Si {dx + F dzf + ^2 dy^] L A = i)dt (3.4) 

z'^ y (1 — ^) J 

where Hi 2-, Si^2, ip and $ are seven functions of x and z to be determined using the numerical 
methods introduced in [Tj. Here, ?/+ parametrizes the black hole temperature as T = y^/An and, 
in the homogenous case, is the black hole horizon size measured in units of the AdS radius. 

The result are hairy, rippled, charged black holes. From the asymptotic form of the scalar field 
one reads off the expectation value of the charged condensate in the dual theory. The condensate 



"'^By adapted here we mean that the T — > hmit can be achieved without introducing large gradients on the 
functions Si 2- 
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Figure 2: The mean value of the condensate as a function of temperature for e = 2 and various 
values of the lattice amplitude. From the inner to outer curves: Aq = 0.0, 0.4, 0.8, 1.2, 1.6, 2.0, 2.4. 
The colors agree with Fig. [Ij 

is a function of x, oscillating about a mean value. In Fig. |2]we plot this mean value as a function 
of temperature for various values of the lattice amplitude. In each case, the condensate follows a 
familiar form, rising rapidly as T drops below Tc and then saturating at low temperature. It is 
clear from Fig. |2]that increasing the lattice amplitude increases the low temperature mean value 
of the condensate. In Fig. |3]we show the variation in the condensate by plotting (0(x)) for Aq = 2 
and various temperatures. The variation is roughly 50% of the mean and the condensate remains 
positive always. 

The transition at T = Tc which turns on the scalar condensate corresponds to a continuous 
phase transition. One can show this without explicitly computing the free energy as follows. 
Consider first the homogeneous case with constant chemical potential fi, and compactify x and y. 
The free energy in the grand canonical ensemble is F = E — TS — fiQ, so its variation is 



5F = 5E- T5S - fi5Q - S5T - QSfi 
But the first law says that 6E = T6S + fi6Q, so at fixed /i. 



dF 



-S 



(3.5) 
(3.6) 



Since the branch of solutions with scalar hair joins the branch of solutions without hair at T 
their entropies must agree, and hence dF/dT is continuous. 

This argument generalizes to the case of a position dependent chemical potential. The /iQ term 
in the free energy is replaced by / fi{x)p{x)dxdy and the fi6Q term in the first law is replaced by 
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Figure 3: The condensate as a function of a; for e = 2, Aq = 2, and T/Tc = .2, .7, .9, .99 from top 
down. Note that by T/Tc = .7 the condensate has almost reached its low temperature limit. 



J fj,{x)6p{x)dxdy (see Appendix). Since we are keeping /x(x) fixed as we change the temperature, 



we recover ( 3.6[ ). The two branches of solutions with ripples still join at T = Tc, so the phase 



transition is again continuous. 

Since we have only introduced the lattice in one direction, one could think of our solutions as 
representing a striped superconductor. However we will not pursue that interpretation herej^ 

In the homogeneous case, the zero temperature limit of these hairy black holes is known to 
take the form j2lj 

(3.7) 



ds' 



r^-dr + dxidx') + 



5for2(-logr) 



and $ = 2(— logr) " near r = 0. This metric has a null singularity at r = 0. The scalar field on 
the horizon of our solutions is becoming more homogenous as T — )■ 0, and at low temperatures, 
the entropy scales like S oc T^ "^ independent of the lattice amplitude. However, the coefficient in 



the entropy formula does depend on the lattice amphtude, so it is not clear if (3.7) applies to the 
T = limit of our rippled solutions. 

The solutions we have discussed so far are not the only static, rippled, hairy black hole solutions. 
There are radial excitations where the scalar field has nodes in the radial direction. There are also 
excitations in the direction of the lattice which can be constructed as follows. Given a solution 
$ with period Aq = 27r//co, — $ is clearly also a solution. One can construct a solution in which 
$ changes sign from one region of size Aq to the next, by simply starting with a seed solution in 

^It is not a gravitational dual of the novel striped superconductor introduced in 23; since that required a 
condensate whose average value was zero. 
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Newton's iterative method that changes sign. Given this, one can construct solutions with period 
(n + m)Ao in which $ changes sign on m of the regions. This would clearly result in a smaller 
mean value for the condensate. We expect that these other solutions all have higher free energy 
than the ones we study, in which $ remains positive everywhere. 

4 Conductivity 

To compute the optical conductivity in the direction of our lattice we introduce a perturbation 
with harmonic time dependence and fix the usual boundary condition on SA^: 

6A^^ h JJx,uj)z + 0(z'^) . (4.1) 

This corresponds to adding a homogeneous electric field = Ee~^'^^ on the boundary. This 
perturbation induces perturbations in most metric components as well as other components of the 
vector potential and scalar. Solving these linear equations with suitable boundary conditions allows 
us to read off the current Jx{x,u) which determines the conductivity via a{u,x) = Jx{x,u)/ E. 
In contrast with the ionic lattice of [8], there are now a total of twelve linear functions to be 
determined. In addition to the perturbations present in the ionic case, we also have two more 
perturbations corresponding to the real and imaginary parts of the charged scalar field. We could 
have worked with U (1) gauge invariant variables, such as in [22], however here we follow a standard 
gauge fixing procedure. For the metric and gauge field perturbations we use the de Bonder and 
Lorentz gauges, respectively. Since we impose a homogeneous electric field, we are interested in 
the homogeneous part of the conductivity, <7{u), that can be retrieved from the full conductivity 
via Fourier mode decomposition: 

a{u,x) = a{u)+ J2 cT„Me'"'=°^ (4.2) 

n=— oo 

In this section, we will fix the lattice amplitude to be Aq = 2, which corresponds to a critical 
temperature Tc = .11/2. 

The real and imaginary parts of the conductivity for T/T^ = .71 are shown in Fig. |4j The first 
thing to note is that there is a pole in the imaginary part of the conductivity. This implies (from 
the Kramers-Kronig relation) that there is a delta function at a; = in the real part confirming 
that this is a superconductor. So, unlike the normal phase, the holographic lattice does not remove 
the delta function in the superconducting regim^ The coefficient of the pole is the superfluid 
density. The bump in the real part of a at w/T 20 is a resonance which is analogous to the 
one studied in [8] coming from a quasinormal mode of the black hole. We will return to this later, 
but for now, the most interesting feature of Re (a) is the rise at low frequency. This shows that 
there is a normal component to the conductivity even in the superconducting regime. Thus, our 
holographic superconductor resembles a two fluid model with a normal component as well as a 
superfluid component. 

^This was noticed earlier in a perturbative treatment of the lattice [TBI I25j . 
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Figure 4: The real and imaginary parts of the conductivity for Aq = 2, /cq = 2, and T/Tc = .71. 
4.1 Low frequency region 

Let us now focus on this low frequency region. The real and imaginary parts of the conductivity 
for various temperatures are shown in Fig. [s] There are four curves ranging from T/Tc = 1 down 
to T/Tc = .7. It should be emphasized that the curve with T/Tc = 1 was done with a completely 
different numerical code using a different metric ansatz and solving for one fewer unknown function 
(no scalar). The continuity of the results is a good check on the numerical accuracy. It is clear 
from Fig. [5] that as we lower the temperature, the normal component is decreasing, and the pole 
in the imaginary part is increasing, showing an increase in the superconducting component. 

A closer examination of the normal component of the conductivity when T < Tc shows that it 
behaves very much like the conductivity above the critical temperature. In fact, at low frequency 
both the real and imaginary parts of the conductivity are very well fit by simply adding a pole to 
the Drude formula: 

aiuj) =% h -. — (4.3) 

bj 1 — IbJT 

where ps is the superfiuid density, pn is the normal fluid density, and r is the relaxation time. 
The three parameters ps, pn and r are temperature dependent but frequency independent. On the 
right hand side of Fig. [5] one can clearly see the Drude behavior in the bottom curve which is still 
in the normal phase. As the temperature is lowered, the pole grows rapidly and soon swamps the 
Drude behavior of the normal component. 

The temperature dependence of ps and p„ are shown in Fig. [6j On the left one sees that ps 
rises rapidly as T drops below T^, and on the right, one sees that p„ drops rapidly. The red dashed 
line is a fit to 

p„(T) = a + 6 e"^/^ with A = 4.0 . (4.4) 
This is similar to a BCS superconductor in which the normal component consists of thermally 
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Figure 5: The low frequency part of the conductivity for T/Tc = 1.0 (blue circles), .97 (red squares), 
.86 (yellow diamonds), and .70 (green triangles). The vertical red line in Re((T) denotes the zero 
frequency delta function. 
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Figure 6: The density of the superfiuid and normal components as a function of temperature. 



extracted from the fit to (4.3). The dashed red line on the right is a fit to (4.4) 
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excited quasiparticles with a gap A, but in BCS theory the gap is smaller, A ^ 1.7 Tc- Observations 
on the cuprates indeed show a gap of order A = 4.0 Tc [26], showing that they are not weakly 
coupled superconductors. In the early work on holographic superconductors, an attempt was 
made to measure this gap. Since the system was homogeneous, there was no Drude peak so people 
looked at \\m^^^Q'Re{a) . This was found to behave like e~^''l'^ over a range of low temperatures. 
In the probe limit (where the spacetime metric is fixed), one found a value close to what we find 
here, A^. = 4.2 Tc [I0l[27], but when backreaction was included, A^- became smaller and depended 
on the charge of the scalar field [H]. With the lattice, we have a much better way to measure 
the superconducting gap and its intriguing that we get a realistic value even for a charge two 
condensate 13 



Another key difference from BCS is the presence of the constant a in (4.4). If we believe this 
extrapolation then p„ remains nonzero even at zero temperature indicating uncondensed spectral 
weight]^ Observations on the cuprates indeed show uncondensed spectral weight at T = [29] . 
Of course the cuprates are d-wave superconductors in which the gap has nodes on the Fermi 
surface which probably contribute to this effect. It is surprising to see p„ 7^ at T = in an 
s-wave superconductor, as we have here. In quasiparticle language, it suggests that the normal 
component consists of two constituents, one of which is gapped and the other remains gapless. Of 
course quasiparticles may not be the right language for the strongly correlated dual system. (A 
precursor to this was seen in the homogeneous case, where it was found that limi^_>o R-e(cr) was 
exponentially small but not zero at T = [24j.) 

To check the extrapolation of p„, one would like to compute pn at lower temperatures. This is 
very difficult because the relaxation time r rises rapidly as the temperature is lowered (see Fig. [?]). 
To see the Drude peak, one needs to probe frequencies cur ~ 1. Since r is rising so rapidly, the 
Drude peak is squeezed to very small frequencies which become difficult to resolve. The dashed 
line in Fig. [7] is a fit to 

r = roe^/'^ with A = 4.3 Tc. (4.5) 

It follows that the scattering rate, 1/r, drops rapidly as T is reduced below Tc. This is another 
observed feature of the cuprates [30]. One often distinguishes clean and dirty limits of supercon- 
ductors by comparing the scattering rate to the gap. It is easy to check that our holographic 
superconductor is in the clean limit: 1/r ^ 2A for all T < T^. 

It is natural to ask if p„ is related to the charge inside the black hole. This appears not to be 
the case. As we lower the temperature, the charge inside the horizon is reduced, and more of the 
charge is carried by the scalar field outside the horizon. However, in the limit T — )■ 0, the charge 



inside goes to zero as the horizon area vanishes, in contrast to (4.4). 



^Curiously, in the probe limit at low temperature, Re(cr) is strongly suppressed for a; < 8Tc, even if one 

changes A^. by changing the mass of the scalar field and spacetime dimension [?7] . 

^In a very different holographic realization of a two fluid model, it was found that p„ remained nonzero at T = 
when the charge on the scalar field was small, but vanished when the scalar charge became larger [25] . 
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Figure 7: The relaxation time as a function of temperature, extracted from the fit to (4.3). The 



dashed curve is a fit to (4.5). 



4.2 Power law 

At shghtly larger frequency, the absolute value of the conductivity of the normal component follows 



the same power law (1.1) as was found in the normal phase above the critical temperature. This 



is shown in Fig. [8] where we have subtracted the pole in the imaginary part of a coming from 
the superfiuid component, and plotted the absolute value of the result, |cr|, minus the off-set C vs 
frequency on a log-log plot. The four lines correspond to the same temperatures as in Fig. [5j The 
fact that the lines are parallel implies that the exponent is the same. The fact that the lines lie on 
top of each other implies that the coefficient B of the power law is again temperature independent. 
Recall that the lowest curve represents the normal phase. This clearly shows that the power law 
fall-off is completely unaffected by the superconducting phase transition. 

Our results are in complete agreement with observations of the bismuth-based cuprate^ [12] . 
Fig. [ojshows measurements of |o"(w) | (without the pole in the imaginary part) for Bi2Sr2CaCu20s+s, 
commonly called BSCCO. The figure includes eight separate log-log plots showing measurements 
on eight different samples ranging from underdoped with Tc = 67K to the optimally doped with 
Tf. = 96K to overdoped with = 60. Each plot shows a variety of temperatures both above and 
below Tc. All show power-law behavior at intermediate frequencies which does not change as the 
temperature is reduced below Tc. The measured exponent is —2/3 for the optimally doped and 
overdoped samples but increases to about —1/2 in the most underdoped samples. The agreement 
between Fig. |8] and Fig. |9] is striking, and we do not understand why our simple gravity model is 
able to reproduce this key feature of BSCCO so well. Of course one important difference is that 



^We thank van der Marel for bringing this paper to our attention. 
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Figure 8: The normal component of the conductivity on a log-log scale showing power law behavior 
with exponent —2/3. The pole in the imaginary part has been subtracted out. The lines correspond 
to T/Tc = 1.0 (blue circles), .97 (red squares), .86 (yellow diamonds), and .70 (green triangles). 



the experimental measurements do not show an off-set C. Another difference is that the power law 
in BSCCO extends up to a temperature independent frequency (the lower cut-off is temperature 
dependent). In our calculations the power law always extends roughly between 2 < ujt < 8 which 
is exactly the same range as was found in the normal phase T > [1]. Since r is temperature 
dependent, both the upper and lower cut-offs are temperature dependent. 

The phase of the complex conductivity computed from the gravitational dual is roughly constant 
over the range of frequencies where the magnitude follows the power law. However, this phase is 
temperature dependent and varies between 60° and 80°. The data on BSCCO shows a temperature 
independent phase of 60°. This difference is likely connected with the constant off-set in our 
power law. Without the offset, scale invariance, causality and time reversal symmetry require 
a(a;) cx {—iu)"' [2J, so the phase is related to the exponent of the power law. 



4.3 Sum rule 

The Ferrell-Glover-Tinkham (FGT) sum rule states that the reduction in the spectral weight when 
T < Tc is taken up by the superfluid density p^: 

/ duRe[aN{i^) - crs{^^)] = ^Ps (4.6) 

where aN is the conductivity in the normal phase and as is the conductivity in the superconducting 
phase. We now ask whether this is satisfied in our gravitational modeQ The left panel of Fig. 10 
shows Re[(T(a;)] over a wide range of frequency in both the normal phase T = and the supercon- 



ducting phase T = .71 Tc. The right panel of Fig. 10 shows that (2/ti) J^J^ duRe[aNi^^) — crsioj)] 



indeed approaches the superfluid density at large u. At the end of the integration, they differ by 
less than 0.6%. This is not only a confirmation that our holographic superconductor satisfies the 



^See [311 1321 133] for general discussions of sum rules from gravity for translationally invariant backgrounds. 
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Figure 9: A log-log plot of the optical conductivity in eight BSCCO samples ranging from un- 
derdoped (UD) to overdoped (OD). The numbers following UD or OD in the plot labels are the 
critical temperatures. Each plot contains curves at several different temperatures both above and 
below Tc. The power-law fall-off is clearly unaffected by Tc. Plot is taken from [T2] . 
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Figure 10: Left panel: The conductivity in the normal phase T = Tc (upper blue curve) and 
superconducting phase T = .71 Tc (lower red curve). Right panel: The red dashed line is the 
superfluid density obtained from the pole in the imaginary part of the conductivity. The black 
line is (2/7r) J^/^ dojRe[(7N{(jj) — crs{uj)]. This integral approaches the red line at large u showing 



that the FGT sum rule (4.6) is satisfied 



FGT sum rule, but also a strong test of the numerics. The two curves on the left panel of Fig. 10 
were computed using different codes, and the fact that an integral of their difference gives precisely 
the expected answer is confirmation of the accuracy of the numerical results. 



Note that the two curves on the left panel of Fig. 10 differ over a range of frequencies of order fi. 
In conventional superconductors, Re[(T(ci;)] is reduced only over a much smaller range of frequency. 



so it is common to cut off the integral in (4.6) at a convenient low frequency ojq. However, when 



this is done in the cuprates, one finds that the integral underestimates the superfluid density [12] . 
In other words, some of the spectral weight represented by ps is missing from the low frequency 
part of Re[cr(u;)]. One must indeed include frequencies of order the chemical potential to recover 
the sum rule, just as we have seen here [M] . 

4.4 Resonance 

The resonance in the conductivity at u/fi ~ 1.5 can be understood on the gravity side as arising 
from a quasinormal mode oscillation of the black hole, just like the example in [S]. As discussed 
in [35] quasinormal modes correspond to poles in retarded Green's functions. We can determine 
the quasinormal mode frequency by fitting the conductivity to 



G^{uj) 1 a + b{u-uJo) 



lUJ 



UJ — UJq 



(4.7) 

For T = Tc, one finds the resonance is very well fit by this expression with Wo/yU = 1.48 — 0.421 
This quasinormal mode changes very little when the temperature drops to .7 T^. 
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5 Discussion 



We have seen that a simple holographic model of a superconductor reproduces quantitative features 
of BSCCO including an intermediate frequency power law with exponent —2/3 which is unaffected 
by the superconducting phase transition, and a superconducting gap A ^ ATc- In addition, it 
reproduces many qualitative features of the cuprates including a rapidly decreasing scattering rate 
below Tc, uncondensed spectral weight at T = 0, and a superconductivity induced transfer of 
spectral weight involving energies of order the chemical potential. 

We find it remarkable that so much of the phenomenology of the cuprates can be reproduced by 
such a simple gravity model. Since gravity coupled to a Maxwell field and charged scalar are part 
of the low energy limit of string theory with many different compactifications, one can perhaps 
view it as the universal part of the theory. In this respect, its dual description might be called 
the "standard model" of strongly correlated systems. One could start with this and then add 
other fields to reproduce observed fine scale structure of particular materials. Perhaps the main 
drawback of such an approach is that our current model describes an s-wave superconductor. We 
still do not have a satisfactory description of a rf-wave holographic superconductor. See [36] for a 
description of some of the difficulties. 

We have presented results for the optical conductivity in the superconducting region for only 
one value of the lattice amplitude Aq = 2 and one value of the wavenumber ko = 2 (with /i = 1). 
We believe that the —2/3 exponent in the power law is independent of these choices since it is 
unchanged when one enters the superconducting regime, and above Tc it has been shown to be very 
robust. We have done preliminary calculations with = 1 which indicate that the behavior of 
the normal component density Pn(T) is also independent of these choices. In particular, although 
Tc increases, the gap remains A = 4Tc and p„ approaches the same nonzero value as T — )■ 0. The 
behavior of the relaxation time r is qualitatively the same, but there is an important difference. 
Although one again finds r = roe^/"^, the value of A is half what it was for k^ = 2. In other words, 
it appears that 1/r oc e~"'^o/T This is the behavior predicted in [T3] where the scattering rate was 
related to the density-density correlation function at a; = and k = k^. Since Poincare invariance 



appears to be restored in the zero temperature infrared geometry (3.7), low energy physical states 
should have cu ~ so states with a; = and k = ko would be exponentially suppressed. 

The fact that we still see -2/3 power law even in the superconducting phase, makes it clear 
that this power law has nothing to do with the AdS2 x R"^ near horizon geometry of the zero 
temperature Reissner- Nordstrom solution. As we have discussed, with the charged scalar present, 
the zero temperature solution has a singular horizon. 

Homes has discovered a remarkable relation between the low temperature superfiuid density 
Ps, the critical temperature Tc and the DC conductivity just above the critical temperature cr^c 
[37] . He found that for a wide range of high temperature superconductors 

Ps^35adcTc. (5.1) 

This relation appears to be universal. It holds for the conductivity perpendicular to the Cu02 
planes in the cuprates as well as the conductivity in the plane, and also holds for the iron pnictides. 
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Unfortunately, it is easy to see that Homes' law cannot hold for the class of holographic supercon- 
ductors we have discussed here where the critical temperature depends on the lattice amplitude]^ 
The reason is that as — ?• 0, adc diverges since the system becomes translationally invariant. 
However, both ps and Tc approach finite limits given by the original homogeneous holographic 
superconductor. So (5.1) cannot hold. To reproduce (5.1) one would have to fix a large Aq and 
vary Tc by other means, perhaps by adding a double trace perturbation |39j. 
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A The first law for black holes with varying chemical po- 
tential 



In this appendix, we derive the first law for electrically charged black hole solutions to (2.1) 
with nonconstant chemical potential, fi{x,y). For convenience we will imagine that the x and y 
directions are periodically identified so the black holes have finite horizon area and mass. We 
follow the approach of Sudarsky and Wald ^O]- As usual for a diffeomorphism invariant theory, 
the Hamiltonian is a linear combination of the constraints plus surface terms: 

H= /"[AT'^Ca + ArM„(A^')] + surface terms (A.l) 

where iV" is the lapse-shift vector and Ca are the usual Hamiltonian and momentum constraints 
of general relativity. The surface terms are determined by the requirement that the variations of 
H with respect to the canonical variables qij,p^\ Ai, E'^ are well defined. In addition to the usual 
gravitational surface terms, there is an additional term coming from the Gauss law constraint: 



/ {N''A^)EidS' (A.2) 

JdT, 



Given a static, electrically charged black hole, we choose S to be a constant t surface which 
starts at the bifurcation surface on the horizon and ends at infinity. We also choose A^'* = (d/dt)"". 
Now under a perturbation of the black hole 



SH = / 6qijdH/dqij + ■ • • = (A.3) 



^An earlier unsuccessful attempt to find a holographic realization of Homes law was made in 
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since the partial derivatives of the Hamiltonian just yield the time derivatives of the canonical 
variables in the background which vanish since the background is static. On the other hand, 
if the perturbation solves the linearized field equations, it certainly solves the linearized con- 



straints, so from (A.l) 6H reduces to a sum of the variation of the surface terms. The gravi- 



tational surface terms yield the usual SM — {K,/8Tr)6A. Since d/dt = at the bifurcation sur- 



face, the electromagnetic surface term (A. 2) only has a contribution from infinity which is simply 
/ y)p{x, y)dxdy. Since y) is fixed by our boundary conditions, the variation of this surface 
term is / fi{x,y)6p{x,y)dxdy. We thus obtain the first law: 

SM = —SA+ / fi{x,y)Sp{x,y)dxdy (A.4) 
Stt J 

In the homogeneous case with constant fi, this clearly reduces to the familiar form 

SM = ^SA + pSQ . (A.5) 
Stt 

References 

[1] G. T. Horowitz, J. E. Santos, and D. Tong, Optical Conductivity with Holographic Lattices, 



JHEP 1207 (2012) 168, |arXiv: 1204.0519 



[2] D. V. d. Marel, H. J. A. Molegraaf, J. Zaanen, Z. Nussinov, F. Carbone, A. Damascelli, 
H. Eisaki, M. Greven, P. H. Kes, and M. Li, Quantum critical behaviour in a high-Tc 
superconductor, Nature 425 (Sept., 2003) 271-274, [ arXiv: cond-mat/0309172| . 

[3] D. van der Marel, F. Carbone, A. B. Kuzmenko, and E. Giannini, Scaling properties of the 
optical conductivity of Bi-based cuprates. Annals of Physics 321 (July, 2006) 1716-1729, 
[arXiv : cond^m at/060403f| . 

[4] A. El Azrak, R. Nahoum, N. Bontemps, M. Guilloux-Viry, C. Thivet, A. Perrin, S. Labdi, 
Z. Z. Li, and H. Raffy, Infrared properties o/YBa2Cu307 and Bi2Sr2Ca„_iCu„02?i+4 thin 
films, Phys. Rev. B 49 (Apr, 1994) 9846-9856. 

[5] P. W. Anderson, Infrared conductivity of cuprate metals: Detailed fit using luttinger- liquid 
theory, Phys. Rev. B 55 (May, 1997) 11785-11788. 

[6] T. Kato and M. Imada, Thermodynamics and Optical Conductivity of a Dissipative Carrier 
in a Tight Binding Model, Journal of the Physical Society of Japan 67 (Aug., 1998) 2828, 
(arXiv:cond-mat/9711208j . 

[7] M. R. Norman and A. V. Chubukov, High-frequency behavior of the infrared conductivity of 
cuprates, Phys. Rev. B 73 (Apr, 2006) 140501. 

[8] G. T. Horowitz, J. E. Santos, and D. Tong, Further Evidence for Lattice- Induced Scaling, 
JHEP 1211 (2012) 102, |arXiv: 1209 . 1098]. 



19 



[9] D. Vegh, Holography without translational symmetry, arXiv: 1301.0537 



[10] S. A. Hartnoll, C. P. Herzog, and G. T. Horowitz, Building a Holographic Superconductor, 



Phys. Rev. Lett. 101 (2008) 031601, [arXiv : 0803 . 3295 



[11] S. A. Hartnoll, C. P. Herzog, and G. T. Horowitz, Holographic Superconductors, JHEP 0812 
(2008) 015, [arXiv : 0810. 1563] ■ 

[12] J. Hwang, T. Timusk, and G. D. Gu, Doping dependent optical properties of 

Bi2Sr2CaCu2 08+7, Journal of Physics Condensed Matter 19 (Mar., 2007) 125208, 
[arXiv : cond-mat/0607653j . 

[13] S. S. Gubser, Breaking an Abelian gauge symmetry near a black hole horizon, Phys. Rev. 
D78 (2008) 065034, [irXiv: 0801 .2977] . 

[14] S. A. Hartnoll and D. M. Hofman, Locally Critical Resistivities from Umklapp Scattering, 



Phys. Rev. Lett. 108 (2012) 241601, [arXiv: 1201 .3917 



[15] K. Maeda, T. Okamura, and J.-i. Koga, Inhomogeneous charged black hole solutions in 

asymptotically anti-de Sitter spacetime, Phys. Rev. D85 (2012) 066003, [arXiv : 1 107 . 3677] . 

[16] Y. Liu, K. Schalm, Y.-W. Sun, and J. Zaanen, Lattice Potentials and Fermions in 

Holographic non Fermi-Liquids: Hybridizing Local Quantum Criticality, JHEP 1210 (2012) 
036, |arXi v:1205.5227j . 

[17] R. Flauger, E. Pajer, and S. Papanikolaou, A Striped Holographic Superconductor, Phys. Rev. 
D83 (2011) 064009, (arXiv: 1010. 1775| . 

[18] J. A. Hutasoit, G. Siopsis, and J. Therrien, Conductivity of Strongly Coupled Striped 



Superconductor, arXiv : 1208 . 2964, 



[19] M. Headrick, S. Kitchen, and T. Wiseman, A New approach to static numerical relativity, 
and its application to Kaluza-Klein black holes. Class. Quant. Grav. 27 (2010) 035002, 
[arXiv: 0905. 1822] . 

[20] P. Figueras, J. Lucietti, and T. Wiseman, Ricci solitons, Ricci flow, and strongly coupled 
CFT in the Schwarzschild Unruh or Boulware vacua. Class. Quant. Crav. 28 (2011) 215018, 
jarXiv: 1104.4489] . 

[21] S. Ganguli, J. A. Hutasoit, and G. Siopsis, Enhancement of Critical Temperature of a 



Striped Holographic Superconductor, Phys. Rev. D86 (2012) 125005, [arXiv: 12 05.3107|. 

[22] G. T. Horowitz, J. E. Santos, and B. Way, A Holographic Josephson Junction, 
Phys.Rev.Lett. 106 (2011) 221601, [a rXiv: 1101 .332 6|. 

[23] E. Berg, E. Fradkin, S. A. Kivelson, and J. Tranquada, Striped superconductors: How the 
cuprates intertwine spin, charge and superconducting orders, .arXiv : 0901 . 4826, 



20 



[24] G. T. Horowitz and M. M. Roberts, Zero Temperature Limit of Holographic 
Superconductors, JHEP 0911 (2009) 015, [arXiv: 09 08. 3677 1. 

[25] N. lizuka and K. Maeda, Towards the Lattice Effects on the Holographic Superconductor, 
JHEP 1211 (2012) 117, I arXiv: 1207.2943]. 

[26] K. K. Gomes, A. N. Pasupathy, A. Pushp, S. Ono, Y. Ando, and A. Yazdani, Visualizing 
pair formation on the atomic scale in the high-Tc superconductor Bi2Sr2CaCu20s+7, Nature 
447 (May, 2007) 569-572, [arXiv:0706.0214] . 

[27] G. T. Horowitz and M. M. Roberts, Holographic Superconductors with Various Condensates, 
Phys.Rev. D78 (2008) 126008, |arXiv: 0810 . 1077] . 

[28] J. Sonner and B. Withers, A gravity derivation of the Tisza-Landau Model in AdS/CFT, 
Phys.Rev. D82 (2010) 026001, [arXiv: 100 4.2707 . 

[29] J. Orenstein, Optical Conductivity and spatial inhomogeneity in cuprate superconductors, in 
Handbook of high-temperature superconductivity. Theory and experiment. (Schrieffer, John 
Robert and Brooks, James S., ed.). New York, NY: Springer, 2007. 

[30] D. A. Bonn, R. Liang, T. M. Riseman, D. J. Baar, D. C. Morgan, K. Zhang, P. Dosanjh, 
T. L. Duty, A. MacFarlane, G. D. Morris, J. H. Brewer, W. N. Hardy, C. Kalhn, and A. J. 
Berhnsky, Microwave determination of the quasiparticle scattering time in YBa2Cu3 06.95, 
Phys. Rev. B 47 (May, 1993) 11314-11328. 

[31] D. R. Gulotta, G. P. Herzog, and M. Kaminski, Sum Rules from an Extra Dimension, JHEP 
1101 (2011) 148, [arXiv: 1010.4806 . 

[32] W. Witczak-Krempa and S. Sachdev, The quasi-normal modes of quantum criticality, 
Phys.Rev. B86 (2012) 235115, |arXiv: 1210.4166| . 

[33] W. Witczak-Krempa and S. Sachdev, Dispersing quasinormal modes in 2-hl dimensional 
conformal field theories, arXiv : 1302 . 0847. 

[34] A. V. Boris, N. N. Kovaleva, O. V. Dolgov, T. Holden, G. T. Lin, B. Keimer, and 

G. Bernhard, In-Plane Spectral Weight Shift of Charge Carriers in YBa2Gu306.9, Science 
304 (2004), no. 5671 708-710, 

[http : //www . sciencemag . org/content/304/5671/708 . full . pdf . 

[35] D. T. Son and A. O. Starinets, Minkowski space correlators in AdS / CFT correspondence: 
Recipe and applications, JHEP 0209 (2002) 042, |arXiv:hep-th/020505lj . 

[36] F. Benini, G. P. Herzog, R. Rahman, and A. Yarom, Cauge gravity duality for d-wave 
superconductors: prospects and challenges, JHEP 1011 (2010) 137, [arXiv: 1007 . 1981 . 



21 



[37] C. C. Homes, S. V. Dordevic, T. Valla, and M. Strongin, Scaling of the superfluid density in 
high-temperature superconductors, Phys.Rev. B 72 (Oct., 2005) 134517, 
[arXiv: cond-mat70410719j . 

[38] J. Erdmenger, P. Kerner, and S. MuUer, Towards a Holographic Realization of Homes' Law, 
JHEP 1210 (2012) 021, [arXiv: 1 206 . 5305 . 

[39] T. Faulkner, G. T. Horowitz, and M. M. Roberts, Holographic quantum criticality from 
multi-trace deformations, JHEP 1104 (2011) 051, [arXiv: 1008. 1581] . 

[40] D. Sudarsky and R. M. Wald, Extrema of mass, stationarity , and staticity, and solutions to 
the Einstein Yang-Mills equations, Phys.Rev. D46 (1992) 1453-1474. 



22 



